Quantum Instability of Magnetized Stellar Objects 



R. Gonzalez Felipe 1 , H. J. Mosquera Cuesta 2 ' 3 , A. Perez Martinez 4 and H. Perez 
Rojas 4 

1 Departamento de Fi'sica and Centro de Fi'sica Teorica de Particulas, Instituto Superior 
Tecnico, Av. Rovisco Pais, 1049-001 Lisboa, Portugal 

2 Instituto de Cosmologia, Relatividade e Astrofi'sica (ICRA-BR), Centro Brasileiro de 
Pesquisas Fisicas, Rua Dr. Xavier Sigaud 150, CEP 22290-180, Urea, Rio de Janeiro, RJ, 
Brazil 

3 Centro Latino-Americano de Ffsica, Avenida Wenceslau Braz 71, CEP 22290-140 fundos, 
Botafogo, Rio de Janeiro, Brazil 

4 Instituto de Cibernetica, Matematica y Fisica, Calle E No. 309, 10400 La Habana, Cuba 



Abstract The equations of state for degenerate electron and neutron gases are 
studied in the presence of magnetic fields. After including quantum effects to study 
the structural properties of these systems, it is found that some hypermagnetized 
stars can be unstable based on the criterium of stability of pressures. Highly mag- 
netized white dwarfs should collapse producing a supernova type la, whilst su- 
perstrongly magnetized neutron stars cannot stand their own magnetic field and 
must implode, too. A comparison of our results with a set of the available obser- 
vational data of some compact stars is also presented, and the agreement between 
this theory and observations is verified. 

Key words: Dense matter — equation of state-instabilities-magnetic fields- stars: 
neutron — white dwarfs 



1 INTRODUCTION 

Many stellar objects are known to be endowed with large magnetic fields (Dryzek et al. 2002 
Heyl 119991 I2000fl . For instance, white dwarfs with surface magnetic fields whose strengths 
range from 10 5 to 10 9 G have been discovered (see an earlier set of references in Kemp et 
al. HT?7T)1 Putney ITMBI Schmidt & Smith ITMBI Reimers ITSMl and the updated list of Suh & 
Mathews 2000). Moreover, magnetic field strengths of the order of 10 20 G have been suggested 
to exist in the core of neutron stars or pulsars. Recently, in Refs. Chaichian et al. 2000 and Perez 
et al. 2003] electron and neutron gases in a strong magnetic field were considered with the aim 
to study the equation of state of white dwarf stars, neutron stars and its relation to supernovae. 
It was found a pure quantum effect: the appearance of a ferromagnetic configuration in the 
neutron star interior, which is intrinsically related to the presence of the magnetic field. Such 
an effect opens the possibility of a quantum magnetic collapse of the gas under consideration. 
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This effect is related to the density of the star and its magnetic field, and as such, it allows one 
to establish a criterium of stability taking into account these physical properties. 

The presence of a magnetic field drives the loss of the rotational symmetry of the particle 
spectrum, which in turn manifests as an anisotropy in the thermodynamic properties of the 
system. This behavior can be seen through the energy-momentum tensor 7^„, if we recall that 
the external magnetic field TL induces a magnetization M in the medium which is described 
through the relation TC = B — 4irA4(B). Here B is the microscopic magnetic field, which is 
assumed to point along the x% axis in what follows. Starting from the energy-momentum tensor, 
one can derive the expressions for the pressure components, longitudinal (p 3 ) and perpendicular 
to the magnetic field, 

T ± = P± = -n-BM , r 33 = P3 = -n , (1) 

where f2 is the thermodynamic potential. 

In classical electrodynamics (Landau & Lifshitz 1970), i.e., when the spin interactions are 
not taken into account, this anisotropy appears but the results are quite different. Since 

B 2 B 2 
P± = Po + -r- , P3 = Po - -— , (2) 

47T 47T 

where po is the isotropic pressure, one obtains at the classical level 

P± > P3 ■ (3) 

This fact explains the oblateness of some astrophysical objects when one studies them in a 
classical way. Proper examples of this are provided by Shapiro & Teukolsky 1983 in the case 
of magnetic white dwarfs, and by Cardan, Prakash & Lattimer 200 J for the highly magnetized 
neutron stars. Contrary to these well-known effects, the interplay of magnetic and quantum 
effects imply that some stars get a cigar-like shape along the X3 axis, and some of them may 
even collapse. 

The purpose of this paper is to discuss in more detail and to exploit some of the astrophys- 
ical consequences of previous works (Chaichian et al. 120001 Perez et al. |2003| regarding the 
appearance of hydrodynamic instabilities in strongly magnetized electron and neutron gases as 
realizations of the physics taking place inside white dwarfs (WDs) and neutron stars (NSs), 
respectively. In particular, these instabilities on the configuration of any super critically mag- 
netized stellar object are shown to appear due to the action of quantum-mechanical effects, as 
the occupation of the particle Landau ground state, which is driven by the cooperative particle 
spin-magnetic field coupling. These new effects allow us to introduce a new model to show 
that some WDs may become ultramagnetized and may collapse in a kind of SN type la, even 
without reaching the Chandrasekhar mass limit. As a matter of consistency, we test this theo- 
retical argument in a wide and extensively studied set of astrophysical sources and show that 
the conclusions drawn from our analysis are in complete agreement with present observations 
of compact remnant stars, while forbid the existence of some exotic configurations of them. 

The structure of the paper is as follows: In Section[21we present the physics of a magnetized 
electron gas and apply it to the study of the stability conditions of white dwarfs. Based on 
this physics, a new model for SN type la events is presented, which not necessarily depends 
upon the WD to reach its critical mass. Section [21 considers the case of a neutron g as, as a 
model for NSs. We show that the so-called magnetars cannot in principle be formed based 
on the instability criterium of anisotropic pressures, while canonical pulsars would. Finally, in 
Section 01 we present our conclusions. 
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2 MAGNETIZED ELECTRON GAS 

In this section we shall study the thermodynamic properties of the degenerate electron gas 
in very intense magnetic fields and density regimes, which for astrophysical scenarios are the 
more interesting properties to be investigated. We shall find the conditions that may lead to 
the vanishing of the transverse pressure of the electron gas, as a model of a white dwarf, which 
may then undergo a gravitational collapse. 

We start by defining the thermodynamic potential as 1 

Q, = -T\nZ, (4) 

where T is the gas temperature and Z the partition function. In the case of a magnetized 
degenerate electron gas, the sum over Landau levels appear in the electron gas configuration 
and, therefore, the thermodynamical potential can be expressed as 



o 



r-n ; — , , s x + \/x 2 - 1 - 2nB/B, 

Jx 2 - 1 - 2nB B c - (1 + 2nB B c ) x In v , ' 

V V 1 + 2nB/B c 

"(5) 

where we define J7o = em 2 / (4ir 2 h 2 c 2 ) , a n = 2 — 5o n , x = /j, e /m e , B c = m 2 /(ehc), with /i e as the 
chemical potential for electrons and corresponds to the maximum Landau level for a given 
Fermi energy and magnetic field strength. The maximum occupancy of these Landau levels is 
defined as 

Be 
2B 



(x 2 - 11 



(6) 



I[x] denotes the integer part of x. 

The mean density of particles is given by N e = —dil e /dfi e . In the degenerate limit, i.e., at 
zero temperature, one obtains 

B 



N e = N [^-\J2 <W^ 2 ~ 1 - 2nB/B c , (7) 

where No = ml/(ATT 2 h 3 c 3 ). Bearing in mind that the magnetization is defined as M. e = 
—dile/dB , from Eq. (JjjJ one obtains 



o 



r-a ; — , , / x + Jx 2 - 1 - 2nB/B c 

xJx 2 - 1 - 2nB B c - (x 2 + AnB B c ) x In v - _ ' 

I y/1 + 2nB/B c 



where we define A4q = em 2 j (4ir 2 h 2 c 2 ) . We stress that the magnetization, as well as other 
thermodynamic quantities, is a function of B and N e . Hence, one can write the expression of 
Q e in terms of the magnetization as 



n e = -BM e - Op V n In I - ' v I ' - I . (9) 

Since the thermodynamic potential is a measure of the internal pressure in the gas, then 
it becomes especially interesting to discuss the equation of state of the system at this point. It 
is worth to note that the latter receives contributions from the partial pressures of the several 
species involved. By calculating the energy- momentum tensor of this gas (made up of electrically 



' x + \J x 2 - 1 - 2nB/B c 



Notice that we are defining the thermodynamic potential per unit volume. 
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charged particles) one obtains different equations of state for the directions parallel, p$, and 
perpendicular, pj_, to the magnetic field B. They read 

P3 = -^e , (10) 

and 

2e 2 B 2 ^ , fx + Jx 2 - \-2nBlB c 
Pi. = n,*. ^ > n hi , = 

It is readily verified from these expressions that for a positive magnetization the transversal 
pressure exerted by the charged particles in the presence of a magnetic field is smaller than the 
longitudinal one by the amount BA4 e . 

There are some situations where densities and magnetic fields are such that only the first 
Landau level of the particle spectrum is occupied. Eq. shows that when the electrons are 
confined to the Landau ground state, one can write 

P± = -fie - BMe = , (12) 

while higher excited Landau states provide a positive contribution to the pressure. If the domi- 
nant contribution to the pressure comes from the electron gas, then the vanishing of p± means 
that the pressure due to the gravitational force, of order GM 2 /R 4 (R is the radius of the star), 
cannot be compensated, and thus there appears an instability which leads to an abrupt im- 
plosion of the electron gas. It is interesting therefore to find the conditions for the occurrence 
of this confinement to the ground state (n = 0) and, consequently, for the triggering of the 
gravitational collapse of the star. 

In general, from Eq. © one can write 

iV 2 

«^2x 10~ 21 ^f, (13) 

which essentially defines an instability criterium for all magnetized stars in the plane (N e ,B). 
In Fig. n we show the instability region (dotted area) for a magnetized electron gas in that 
plane. Any star for which its structural configuration places it inside the dotted region should 
undergo a transverse collapse due to the vanishing of the pressure transverse (perpendicular) 
to the magnetic field. 




2.1 White Dwarfs 

Recent observations have pointed out that supernovae type la (SNIa), considered as standard 
candles in the redshift range 0.1 < z < 1, seem to provide evidence for our Universe to be 
accelerating (Riess et al. 1998, Perlmutter et al. I1999II . driven by some kind of dark energy: 
a cosmological constant or quintessence field. However, the premise underlying this suggestion 
needs to be firmly established on astrophysical foundations. In order to settle out this issue 
in a conclusive way, the redshift independent nature of SNIa should be demonstrated (Karl et 
al. 2003, Napiwotzki et al. 2003, Hillcbrandt et al. 2003), since only a complete understanding of 
those events can give the required confidence, if any, to the idea that they are indeed standard 
candles. In this paper we introduce a new scenario for the SN type la events based on the physics 
presented above, which indicates that another kind of collapse may produce SNIa outbursts. 

The leading scenario for SNIa suggests that these explosions stem from the complete disrup- 
tion of a white dwarf (WD) induced by accretion of mass from a companion star, which forces 
it to surpass the WD Chandrasekhar mass limit M cr i t ~ 1.445 M (Livio EHOOIl ■ The abrupt 
conversion of nearly 1 M Q of C/O to 56 Ni and the 56 Ni subsequent decay to 56 Fe releases, on 
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Fig. 1 The instability region (dotted area) in the (N e , _B)-plane for a magnetized 
electron gas. A WD star whose configuration lies inside the dotted region collapses 
due to the vanishing of the transverse pressure p± = — £l e — BA4 e . The labelled points 
represented by stars correspond to some of the stable white dwarf configurations given 
in Table □ 



WD 


M/Mq 


N e /N 


B (in Gauss) 


A 


1.28 


1.4 


1.3 x 10" 


B 


0.9 


0.6 


2.3 x 10 6 


C 


1.2 


1.2 


3 x 10 6 


D 


1.31 


1.7 


3.5 x 10 6 


E 


0.76 - 1.0 


0.3-0.9 


> 3 x 10 8 


F 


> 1.0 


0.9 


3.2 x 10* 


G 


1.35 


1.9 


6.6 x 10 s 


H 


? 


? 


~ 10 9 


I 


? 


? 


> io 9 



Table 1 Mass M, electron number density N e and surface magnetic field B for some 
typical magnetic white dwarfs: A = PG 0136 + 251, B = PG 2329 + 267, C = 1RXS 
J0823.6-2525, D = PG 1658 + 441, E = LB 11146B, F = Grw + 70°8247, G = RE 
J0317-858, H = PG 1031 + 234 , I = GD 229. The above data were taken from Suh & 
Mathews I3UUU1 Here iV = m 3 J(4ir 2 h 3 c 3 ) ~ 5.9 x 10 29 cm~ 3 and M Q ~ 1.99 x 10 33 g 
is the solar mass. 



the deserved timescale, an amount of energy just as large as needed to power the observed lu- 
minosity and kinematics of a SNIa. In addition, this scenario provides a consistent explanation 
for the absence of hydrogen in these events. Although this model is well motivated, it has been 
shown to contain several drawbacks that tend to disfavor it (Hillebrandt et al. [2003 ) . 

In an alternative view, recent Hubble Space Telescope observations suggest that a good can- 
didate to SN type la progenitors could be coalescing binary WD- WD dubbed as blue strugglers. 
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These objects are a class of stars that live in globular clusters orbiting in the outskirts of spiral 
galaxies. A particular set of these candidates to SNIa are mergers of unequal mass WD- WD 
binaries that come together because of angular momentum loss via gravitational radiation (GR) 
emission (Karl et al. EZUUl Napiwotzki et al.|.2003 ). The timescale (~ 10 6 yr) for the system to 
coalesce is dictated by GR, so that the donor companion has the right time scale to fill its Roche 
lobe driving then mass transfer onto its heavier orbital partner fNomoto ll998[l . However, there 
exists a debate in the literature whether an explosion does occur under such a dynamics. Some 
authors contend that view and suggest that more likely a "silent" accretion-induced collapse 
(AIC) to form a remnant neutron star must take place (Mochkovitch et al. I1997[) . In overall, 
the issue seems contentious. 



2.2 Quantum-Magnetic Collapse of a WD as Model of SNIa 

In the lines of the theory presented earlier (Perez et al. l2003[l . we suggest here a quite different 
mechanism in which, for the case of young, rapidly rotating, highly magnetized and massive 
C/O WDs, a direct collapse of a single WD should indeed occur. The overall process is driven 
by both the quantum and magnetic effects described above. 

Rapidly rotating white dwarfs have been studied theoretically (see Shapiro & 
Teukolskv ll983l and the extended list of references therein). Despite earlier thoughts dismissing 
that possibility, based on the fact that the by-the-time observed H lines in most WD spectra 
exhibited narrow cores, recent observations, however, confirm that such objects do exist (Karl 
et al. 12001(1 among the more than 2400 WDs cataloged up today. The shortest rotation period 
that has been clearly measured is the 33 s of AE Aqr, while a good candidate for a slightly 
faster spin is WZ Sge with a period of 28 s. Nonetheless, there appears to exist a faster class of 
WDs spinning with periods of P cr u ~ 27 s, dubbed as DQ Herculis, which, on the other hand, 
are magnetized WD accreting matter in binaries christened cataclysmic variables. 

For non-magnetic white dwarfs, meanwhile, the theoretical view suggests that there exists 
a critical rotation period (or frequency) of uniformly spinning WDs for which mass shedding 
at the star equator takes place. The instability sets in for a Keplerian rotation frequency 

n K = 0.67(7rGp) 1/2 , (14) 

where p is the WD mean density. That relation implies a rotation period Pk ~ 17 s, for stars 
with mean density in the range 10 4 g/cm 3 < p < 10 7 g/cm 3 . However, for WD stars rotating 
differentially at the relevant temperature the threshold period can be much smaller than the 
above limit. This is a crucial issue that we will take advantage of in developing our model. 
According to Shapiro & Teukolsky[t983 in a differentially rotating WD the ratio between the 
frequency at the center and at its surface is 

^equator + 1 /-i r-\ 

O ~K> ( 15 > 

& ^center 

for almost all the equilibrium configurations with masses in the range 0.4 — 0.9 M©. Massive 
WD models above 1.4 Mq possess surface velocities around 3000 — 7000 km/s, which suggests 
periods around 1 s. 

A WD in a binary system may accrete material from its orbital partner in such a way that 
its angular momentum can rise significantly, and it may rotate differentially. As the material 
accreted transfers energy and angular momentum, the primary WD star will be spun-up over 
the time scale 

. , J 2MR 2 fl 1 fi / M x 1/2 
AT ~ — = r- — = 4 x 10 6 yr. 
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(16) 



where J is the angular momentum accreted from the companion, J - the applied torque, P - the 
rotation period, At - the accretion rate, and R is the WD radius in units of the solar radius Rq. 
Notice that a higher accretion rate, as in unequal WD- WD binaries, may reduce the spin-up 
time scale. Therefore, an essential ingredient of this model is the existence of a WD primary 
in a binary system which is old enough to having had time to accrete up to the point of being 
spun-up to its break-up velocity at the equator. 

Once we have the star differentially rotating near its break-up period, the conditions 
are reached for the a — fl dynamo mechanism (Duncan & Thompson 1992; Thompson & 
Duncan 19931 119951 l*HK)6) to start to amplify the WD initial magnetic field. The maximum 
magnetic field achievable through this process can be as strong as (the equipartition of pressure 
and kinetic energy) 

& Mto^pC _> B . „ 3 x 10 i4 f —\ 1 G , (17) 



8tt " (87t/3)i4 D c Vis 

with B s c being the supercritical B-field and Pi the shortest period reached upon undergoing the 
a — fl dynamo phase. This B s c should be generated as the differential rotation is smoothed by 
growing magnetic stresses inside the WD star. For the typical WD mass, radius and rotation 
frequency given above, the associated spin energy in the process reads 

Em = ~M WD f4 D i4 D ~ 10 51 (J^J erg , (18) 

which is similar to the characteristic energy of most SNIa events. This energy plus the gravita- 
tional binding energy, which is of the same order of magnitude as the above, would be released 
during the implosion of the star. The computed amplified magnetic field 2 plus the matter den- 
sity of a typical WD will put the star in the instability region showed in Fig ^ according to 
Eq. H12(l . where it must undergo a gravitational collapse. Thus, the higher the WD initial in- 
field, the shorter the break-up spin-up time scale. This process may take place well in advance 
the star reaches the Chandrasckhar critical mass of 1.445 M@. 

Since the time scale estimated above is shorter than the age of our universe, certainly old 
WD stars may accrete enough material so as to rotate at its break-up velocity. At that point 
its initial _B-field might be enlarged via the dynamo mechanism, which in turn may drive the 
star to explode as a consequence of the quantum-magnetic collapse it undergoes. Those systems 
may be triggering the observed SNIa events. 

Whenever the physical properties of the WD star (see Table ^) places it above the insta- 
bility region depicted in Fig. ^ the WD must implode triggering a supernova-like event with 
the characteristics of a SNIa. The WD dynamical timescale and the abrupt combustion of its 
constituent material (electron capture into protons, and rapid process nucleosynthesis to 56 Ni) 
is to produce the energetics as required on the observational basis. The main reason for this 
claim is that both characteristics, i.e. the WD rapid spinning and age, can amplify (via the 
a — fi dynamo mechanism) the remnant magnetic field (~ 10 s G) up to a critical value to 
which it may be unstable and thus collapses. In fact, there is strong evidence for young and 
highly magnetized WDs, as for instance the object RE J0317-858 listed in Tabled The process 
we describe is alike to the accretion induced collapse quoted above, but originated from the 
cooperative action of two different physical effects. At the end, the SNIa explosion must be 



Notice that the final _B-field is almost independent of the WD initial B-field. 
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accompanied by a gamma-ray burst of energy ~ 10 51 erg. This is an interesting signature to 
look for. A further testable prediction of this picture could be to find no remnant at all after 
a SNIa! Such a result could be interpreted as if the WD had collapsed into a black hole or a 
black string, driven by the combination of the quantum and magnetic effects discussed here, 
and purported in Perez et al. 120031 

Were this model, i.e. the collapse of heavy and strongly magnetized WD, realized in nature, 
it would significantly help to clarify whether SNIa could indeed be considered as legitimate 
standard candles, since for the collapse to take place in this case the WD mass must be in 
the range 0.9 — 1.35 M©, for which the total binding energy released is practically the same 
~ 10 51 erg. Of course, some other observational constraints play some role, as for instance the 
statistics of massive WD- WD binaries, which is relatively small (Karl et al. 120031 Napiwotzki 
et al-EHSJ. 

In brief, several white dwarfs endowed with strong magnetic fields have been observed 
and studied (Kemp et al. 119701 Suh & Mathews EQH01 see Table [TJ. From the analysis of the 
observational data (Barstow et al. [TSHSI Greenstein et al. I19820 one can see that our model is 
in close agreement with those observations. These stars are stable objects. On the other hand, 
if those WD listed in Table ^were given an amplified supercritical (B%) magnetic field, then 
they must implode. 



Most of the observed neutron stars are pulsars, i.e., fast rotating neutron stars with strong 
magnetic fields (Lorimer 1999J. They consist mainly of neutron matter with a high central 
density. As in the case of WDs, a similar study can be performed for neutron stars which are 
properly described by a degenerate neutron gas (Shapiro & Teukolsky[l983 ). In this case, it is 
straightforward to calculate the thermodynamic quantities starting from the neutron density 
N n . Assuming that the uniform magnetic field B lies along the 23-axis, and using the Dirac 
equation for neutral particles with anomalous magnetic moment (Bagrov & Gitman 1990) 
propagating through this field, we get the spectrum as 



where p^,p± are the momentum components along and perpendicular to B, respectively; m n 
is the neutron mass and q n ~ — 1.91/zjy is the neutron magnetic moment (/in — e/(2m p ) is 
the nuclear magneton), and 77 = ±1 are the 03 eigenvalues corresponding to the two possible 
orientations (parallel or antiparallel) of the neutron magnetic moment with respect to the 
magnetic field. 

The thermodynamic properties of a degenerate neutron gas are easily obtained following 
the same procedure as in Section II and it was the scope of Ref. Perez et al. |2P03> which we 
revisit now to discuss their astrophysical implications. 

Firstly, the neutron density N n , thermodynamical potential f2 n , and magnetization M n 
read, respectively, 



3 MAGNETIZED NEUTRON GAS 




(19) 




(20) 



77=1,-1 . 




+ 



(1 + r]y)(5i]y - 3)xf v (1 + rjy) 3 (3 - ijy) 



24 24 



rjyx 3 
6 



s 



,(21) 
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M n = -Mo ? ? 

n=i-i 



(l-2r,y)xf n (1 + 77^)2(1 - 77^/2) 



(22) 

Here we define the numerical constants: No = m n / (4ir 2 h 3 c 3 ) ~ 2.73 x 10 39 cm -3 , fio = Nom n ~ 
4.11 x 10 36 erg cm" 3 and Mo 



N q n ~ 2.63 x 10 16 erg/(Gcni ). We also set 



with £>„ 



m n /q r , 



1.56 x 10 G. and introduced the notations: 



/„ = v^ 2 - (1 + W) 2 



7T 

S " = 2 



1 + 173/ 



hi 



X + fr, 

1 + rjy 



(23) 



(24) 



In the limit B = 0, Eqs. H2()J) and (|21|l reproduce the usual density and thermodynamic 
potential of a relativistic Fermi gas at zero temperature (Fradkin 1967). We see that the mag- 
netization M n is a nonlinear function of B, and since M n > the magnetic response of the 
neutron gas is ferromagnetic (Perez et al. 12003) ) . 

As we mentioned in the introduction, the key ingredient brought by the anisotropy of 
pressures is precisely that the magnetization is positive. In fact, one finds that the condition 
P3 > P± holds for any study of a highly magnetized Fermi gas. In particular, we remark that 
for a neutron gas with anomalous magnetic moment this condition is fulfilled. 

We also notice that if x = 1 + rjy, then f v = s v = L n = 0. For x < 1 + y the thermody- 
namic potential f2„, and consequently, all the thermodynamic quantities in the system become 
complex. The condition x = 1 + y defines a curve in the (N n , B)-plane described by the relation 



arcsm 

2 



1 



y 



(25) 



which delimits the region where the pressure becomes complex. Since ?/ <gc 1 for the magnetic 
field strengths under consideration (B s ~ 10 15 G in the NS surface), we can expand Eq. (I25|l 
around y = to obtain the approximate expression 



N n ~ -N y A ". (26) 

It is also straightforward to find the curve which corresponds to the vanishing of p±. Writing 
x = 1 + (l + a)y, with a > in Eqs. (|20|l - l|22|) . and expanding around y — we obtain to leading 
order in y : 

2^2 . 



N n 



-No y 



3/2 



(a + 2) 



3/2 , „3/2 



P± 



2\/2 y 



5/2 



2a 3 + 7a 2 + 4a + a 3/2 (5 + 2a)V 2 + a - 



(27) 



(28) 



The solution of the equation p± = is given by a 
value into Eq. (|27|l wc find 



— 3 



y/E- 1 ~ 0.34. Substituting this 



2\/2 



No y 



3/2 



V5 + 1 



3/2 



'-Vb - 1 



3/2' 



3.57iVo y 



3/2 



(29) 



In Table II we present some results for neutron stars obtained by using other equations 
of state (EOS), and taking into account that B-fields in the core may fulfill the dipolar law, 
i.e., that B obeys the relation B = B s (i? 3 /r 3 ), at least for the dipole component; here R s is 
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NS 


EOS 


M/Mq 


N n /N 


B (in Gauss) 


A 


Pol2 


2.31 


0.019 


5.1 x 10 2U 


B 


Pol2 


2.37 


0.11 


1.48 x 10 20 


C 


BJI 


1.96 


0.18 


9.1 x 10 20 


Y) 


BJI 


2.12 


0.25 


9 86 x 1 21 


E 


PandN 


1.72 


0.26 


5.8 x 10 20 


F 


PandN 


1.86 


0.35 


1.6 x 10 21 


G 


Akmal 


2.22 


0.21 


6.4 x 10 20 


H 


Akmal 


2.53 


0.19 


3.2 x 10 21 


I 


PCL 


2.09 


0.054 


8.7 x 10 20 


J 


PCL 


1.91 


0.25 


1.6 x 10 21 


K 


PCLhyp 


2.02 


0.30 


6.9 x 10 20 


L 


PCLhyp 


1.76 


0.04 


1.38 x 10 21 



Table 2 Mass M, neutron number density N n and magnetic field B for some neutron 
star equilibrium configurations computed in Ref. l200Tl using different equations of state 
(EOS). Here N Q = m 3 /(4tt 2 /i 3 c 3 ) ~ 2.73 x 10 39 cm" 3 . 



the star (surface) radius and r is any inner radius of the star. Our results were obtained by 
using the data from Cardall, Prakash & Lattimer l2()()ll where the authors performed "realistic" 
numerical simulations of strongly magnetized neutron stars. We identify B s = B po i e , as of 
usage in Cardall, Prakash & Lattimer 2001 Under these conditions, and based on the proposed 
criterium of stability, we can conclude that those neutron stars would be unstable objects. 

In Fig. |2 we plot the curve defined by Eq. (|29|l in the (N n , B)-plane. For values of x > 
1 + (1 + ao)y the gas pressure is positive and the neutron gas is stable. One can verify, on the 
other hand, that if x < 1 + (1 + &o)y the pressure becomes negative or vanishes, thus leading 
to the transverse collapse of the star. 

3.1 Magnet ars 

It has been recently suggested that there might be stellar objects, known as soft gamma re- 
peaters (SGR), whose magnetic field inferred from the spindown of their believed associated 
pulsars 3 seems to exceed the critical value: B c — m 2 e /{ehc) ~ 4.41 x 10 13 G. In particular, it 
is claimed that these are newly born neutron stars (Mazets et al. E)7lK Murakami et al. 1984 
Kouveliotou et al. I1999[l . also called magnetars (Duncan & Thompson 1992, Thompson & 
Duncan HT9931 ITM51 HD96 1 . with surface magnetic fields as strong as 10 15 G. According to 
Duncan & Thompson: 1992;, Thompson & Duncan 119931 IT§§51 H 996 a neutron star with period 
Pj ~ 1 ms can support an efficient a — f2 dynamo. Since neutron stars are formed with significant 
differential rotation, the associated energy is E^i ~ 10 52 (Pi/l ms) -2 erg. Thus magnetic fields 
as strong as 3 x 10 17 (P;/1 ms) -1 G can be generated as the differential rotation is smoothed 
by growing magnetic stresses inside the star. After the available energy is released in the outer- 
most parts of the star, vigorous convection continues to generate much stronger magnetic fields 
than any previous phase of stellar convection, producing fields larger than the one required 
to balance the gravitational binding energy density (Duncan & Thompson H9521 Thompson & 
Duncan 1993, 19951 1199!))) . Thus their result suggests that the a — fl dynamo operating in 1 ms 
neutron stars might generate surface dipole fields largely stronger than 10 13 G. 



3 See also Guesta l^OOOI for a different approach. 
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Fig. 2 The instability region in the (N n , B)-plane for a magnetized neutron gas, as 
a model of a pulsar. A star whose configuration lies above the solid curve defined 
by Eq. I)29|l should collapse due to the vanishing of the transverse pressure p± = 
—Cl n — BM. n . The labelled points represented by stars correspond to the neutron 
star configurations listed in Tableland computed using different EOS. The core of 
any star being trapped in the instability region must implode or collapse as in the 
quark-nova model (Ouyed et al. 2002). 



For this class of compact magnetized stars with B s ~ 3 x 10 17 (Pi/l ms) _1 G, our description 
above predicts that such objects should implode to a magnetically (B < B c ) stable configura- 
tion, which is not a magnetar! This result, previously obtained in Perez et al. 120031 was later con- 
firmed by Khalilov (Khalilov 2002 ) and Chakrabarty's group in a very interesting series of papers 
(Ghosh & ChakrabartylHIl Ghosh et al-EUl Mandal & Chakrabartv l^OO^al l2002b1 12OT3|) . In 
particular, Khalilov in his stability analysis of a degenerate neutron gas in chemical equilibrium 
with a background of electrons and protons, also included the contribution from the anoma- 
lous magnetic moment of the fermions composing the star, as done in Ref. Perez et al. 2003 
We notice, however, that in Ref. Khalilov 2002 a longitudinal collapse is found instead of a 
transversal one. A detailed discussion of the approach followed in the latter work, which seems 
to us inappropriate, will be presented elsewhere (Perez et al.[5005). 

As commented above (cf. Fig. |2J some neutron stars described by other EOS producing 
B s ~ 3 x 10 17 G fall in the region of instability defined by our theory. The attentive reader 
should note that this magnetic field strength is identical to the one invoked by Kondratyev 
(Kondratyev 2002 ) in an attempt to explain the statistical properties of magnetars. By relating 
intervals of intense activity with sharp, step-like changes of magnetization, due to structural 
inhomogeneities in the crust and release of its stored magnetic energy, Kondratyev's letter 
claims that this could be the mechanism triggering SGRs activity. The demagnetization takes 
place as random jumps associated with magnetic avalanches and sharp energy injection into the 
magnetosphere develop. Because the nucleons populate discrete energy levels in such a matter, 
for field strengths B ~ 10 16-17 G (and density near p n ucicar), where the energy levels crossover, 
the nuclei structure changes abruptly to a stepwise field dependence of the nucleus magnetic 
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moment on the nucleon magneton. According to this author, the decay of this configuration 
would then drive the SGRs bursting phases. 

However, for a crustal multipolar field configuration one can expect, at the stellar interior, 
a field strength much higher than B ~ 10 17 G, at least for the dipole component. This is an 
idea the astrophysical community is well acquainted with. If this is the case, then our analysis 
jumps into scene showing (see Fig. |2J) that much earlier in the neutron star life it should 
have collapsed. Notice further that Woods et al. 120021 after studying large torque variations in 
two SGRs, concluded that within the context of the magnetar model seismic activity cannot 
account for both the bursts and long-term torque changes unless the seismically active regions 
are decoupled from one another. The idea is that since the observed changes in spindown rate 
do not correlate with burst activity then the physical mechanisms behind either phenomena 
are more likely unrelated. Whether the magnetar model is able to cope with this observational 
requirement is not apparent. 

4 CONCLUDING REMARKS 

If a hypermagnetized neutron star could somehow be formed in a supernova explosion, the 
abrupt amplification of its magnetic field will drive it into collapse. Since magnetic flux is 
dissipated during the implosion via a process analogous to the Sun's coronal mass ejections (as 
in Malheiro et al. I2U(J4|) . the most likely outcome of such a collapse may give as a remnant a 
strange star of canonical magnetic field, black hole or even an exotic black string. Once again, 
the violent "consumption" of the neutron star nuclear material may trigger the emission of a 
gamma-ray burst of overall energy ~ 10 52 erg, inasmuch as in the quark-nova model (Ouyed et 

al. ESSll- 
In the case of neutron stars such as the suggested magnetars, the anisotropy of pressures 
must naturally develop, and the condition for the collapse in the direction perpendicular to the 
dipole magnetic field could be satisfied (Fig. |2J) for the typical values of density and magnetic 
field strength routinely quoted for these stars. 

To summarize, we have presented a consistent theory to discuss the stability of compact 
remnant stars whose structure is dictated by a combination of quantum and magnetic effects. 
We have shown that, in general, the theory agrees quite well with current observational data for 
magnetic WDs (Shapiro & Teukolsky 1983}, canonical pulsars (Lorimcr 1999), and the realistic 
relativistic models of Cardall, Prakash & Lattimer 2001. A major outcome is the fact that some 
special configurations of highly magnetized white dwarfs could collapse triggering explosions 
similar to a SNIa. Besides, some neutron stars endowed with superstrong magnetic fields are 
shown to be naturally unstable, and therefore should collapse, for their quoted surface magnetic 
fields. This would drive a powerful supernova explosion followed by a gamma-ray burst, too. 
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